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We consider a self-interacting scalar field in a de Sitter background and deal with the associated 
infrared divergences in a purely diagrammatic way using the in-in formalism. In the particular 
case of a large N 0(N) invariant scalar field theory with quartic self- interactions we recover the 
result that the connected four-point correlation function, which is a signal of non-Gaussianity, is 
non-perturbatively enhanced with respect to its tree-level value. 



00 ■ I. INTRODUCTION 

o ■ 

\ The detection of non-Gaussianity (NG) [l[ in the cosmological perturbations generated during inflation [5| has become 
J ^ ' one of the primary targets of forthcoming experiments measuring the properties of the Cosmic Microwave Background 
O (CMB) and matter fluctuations. NG is originated by the self-interactions of the flelds involved in the inflationary 
, dynamics plus those induced by gravity. The effect of loop corrections on the physical observables generated during 

^ ■ a de Sitter epoch of exactly exponential exp ansion have attracted much attention in the past 0, 4, 5, 6, 7, d, and 

more recently [13, IHI, [13, lH, El, E [H, H E E S lU, HI , especially as far as the resummation of infra-red (IR) 
divergences is concerned. IR divergences appear because of the cumulative effects of the superhorizon perturbations 
r— 1[ which are continuously generated during the de Sitter stage. Different approaches have been put forward to deal 
^ with them. The 2PI (Two-Particle-Irreducible) formalism has been adopted in Refs. [1, 01; while a stochastic fleld 
1^' theory method has originally been used in Ref. [3] where the underlying idea is that the IR part of a scalar fleld 
wp. niciy be considered as a classical space-dependent stochastic fleld satisfying a local Langevin equation. The stochastic 
noise terms arise from the quantum fluctuations which becomes classical at horizon crossing and then contribute 
to the background. An hybrid method, combining the stochastic approach and the out-of-equiHbrium fleld theory 
CN ■ techniques of the in-in formalism to solve the gap equation descriving the time-dependent evolution of the two- 
<^ ' point correlator, has been recently used in Ref. [2l[ for a self-interacting 0{N) model in the limit of large TV and it 
' was shown that the connected four-point correlator, the so-called trispectrum, is non-perturbatively enhanced with 
[ respect to its tree-level value. 

. In this paper we adopt a purely diagrammatic approach based on the in-in formalism to analyze the IR diver- 
• ' gences. We restrict ourselves to a quartic self- interacting scalar fleld in a de Sitter background and show that the 
,: IR resummation of a certain class of diagrams occurs. Generali zing the computation to N scalar flelds subject to an 
QQ ' 0{N) symmetry with large iV, we recover the results of Ref. [2l|. In this sense, our results should be considered 
complementary to alternative approaches, e.g. the stochatic approach. 
• The paper is organized as follows. In Section [ll] we summarize the in-in formalism and the Feynmann rules needed 
', to calculate the higher order corrections for the scalar fleld correlators. In Section Hill we explicitly perform the loop 
^ calculations. In Section IV, we discuss the particular case of 0{N) symmetric model and the trispectrum. In Section 
5h |V]we summarize and conclude our work. 

a 

_ _ I 

II. SELF-INTERACTING SCALAR FIELD IN DE SITTER BACKGROUND 

A. Closed Time Path Formalism 

We use the in-in formalism, also dubbed Closed Time Path (CTP) formalism, to calculate the correlation functions 
of a scalar fleld in a de Sitter background. Following Schwing er [IS HI, we introduce the two external sources J+(.t) 
and J~ (x) and consider the quantity 

Z[J+,J-] = j-{0^\Q+)j+. (1) 
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The vacuum evolves indipendently under two sources J+ and J . We may rewrite the latter as 



Z[J+,J-] 



Texp 



Texp 

t* 



dt / d^xJ {x)(j)H{x) 



dt I d^xJ'^{x)(j)H{x) 



(2) 



where T denotes antitemporal order. Here {(p) is an element of a complete, orthonormal set of common eigenvectors 
of the field operators at some late time t* , 



From the definitions H]) and one can obtain the following relations [2^: 

Z[J,J] = 1, z[j+,j-]= {Z[J-,J+]) 

and 

9"+"Z[J+, J-] 



(3) 
(4) 



aj-(xi) . . . dJ-ix^)dJ+iyi) . . . dJ+{yn) 

{0-\f[(j>H{xi) . . . q^HMmniyi) ■ ■ ■ (l)H{yn)]\0- 



(5) 



The expectations value can be be obtained by variation of the sources J+ and J . In particular for a time-dependent 
Hamiltonian system H{t) that starts in a state \in) at time ii„, we can write the expectation value as: 



(Qit)) 



Texp i / dt' H{t') 



T exp 



dt' H{t') 



(6) 



Now we move to a curved space, namely to a de Sitter background. We write the Lagrangian density for a scalar field 
with potential V(4>) as 



(7) 



where the metric has signature — V ++, ^ is the conformal parameter and is the counterterm. Choosing m = 
and ^ = we select a massless minimally coupled scalar field. The generating functional becomes 



Z[J+,J_,p(iin)] = / I?0+I?^r(0+|p(ii„)|0i-) 

' X>(j,+V(l)-e -^''n ^ <i''x{c[<t>+]-c[^-]+J+<p++J^'l>- ) 

The path integral on the second fine can be written in short-hand notation as 

V(j) exp "i J dt' J d^x (£[</>] + J0) 



(8) 



(9) 



where C is the so-called Schwinger-Keldysh contour which runs from tm to t and back. The field and source J are 
split up in J+ on the first part of this contour, and J_ on the second part, with the condition = (l)^{t). 

To calculate perturbatively correlation functions we need to have the free two-point functions. There are four possible 
time orderings and, using eq. ^ one obtains: 



G-+{x,y) = i{<t>{x)4>{y)Y^) , 
G+-{x,y) =z(0(y)(/.(x))("). 



(10) 
(11) 
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""^ Re(T) 
FIG. 1: An example of Keldysh Contour C 

and 

G++{x,y) = i(T(/)(x)0(2/))(°) (12) 

= ^'(a^o - yQ)G~~^{x,y) + 0{yo ~ xo)G^^{x,y), 
G—{x,y) = j(T0(x)0(2/))W (13) 

= 0{xo - y(i)G^'{x,y) + Oijjo - xo)G^^{x,y), 

where by the superscript (0) we denote the free field correlation functions. They obey the important identity 

G++{x, y) + G—{x, y) = G-+{x, y) + G+- (x, y), (14) 
and they can be put together in a matrix: 

Note that the two point functions depend on the initial conditions via the dependence on p{ti) of the generating func- 
tional eq. It is useful to transform the and 4>~ fields to a different basis, which is a variation of the Keldysh basis: 



c \ _ /(0+ + r)/2 



R I- > (16) 



with 



R.(fl/f|, ,17) 



The free two point functions in this basis can easily be obtained by the transformation 

G.[x, y) = RG(., ,)R- = ( ^^^g'j) ^-^o^' ) , (18) 

The "Gaa" propagator in the matrix l|18p (the element (2,2) of GK{x,y)) is identically zero due to eq. lfT4|) . as can 
be seen by performing directly the product. Finally the Gr and Ga two point functions are called the retarded and 
advanced propagators and GA{x,y) — Gii{y,x). 

B. Feynmann Rules for 

We choose the potential to be V{(j)) = j\4''^- The Lagrangian density becomes: 



-9 (^gf^'^ld,,4.d,^+^cf>^ 



We perform now the field transformation as in eq. lfT6|) . obtaining 



(19) 



We notice that the theory has two vertices. From now on we will utilize the conformal time t, defined as r = 
- dt'/a{t'). Notejhat, as a function of r, the scale factor is a(r) = —(Ht) ^. The free two-point functions in the 
late time limit are [l2i |: 



G'^'ik,n,T,) = 2fc3' (20) 



'1' '2; - 2fc3' 

GRik,Ti,T2) = eiri - T2)^iTi' ~ T2^) (21) 
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and Gfik,Ti,T2) = Gfik, T2,Ti). The two point functions depend only on the length of the spatial momentum 

Following [121 , we represent the (j)c field with a full line and the 0a field with a dashed line and so we can write the 
Feynman rules for the two-point functions as 



Gg'(fe,ri,r2) 



-iGg'(fc,ri,r2) = -iG^l\k,T2,ri) . 

We have two different vertices. One contains three powers of (fic and one of (fiA so we draw it with three full lines 
and one dashed Hue. The other instead contains one power of ipc and three of 0a, hence a vertex with three dashed 
lines and one full line. Since we are in a de Sitter background, y/—g — a'^i'T) and the vertices become: 





-za4(T)A0c'0A -ia%T)^cj,ccljl 
When a two point function is attached to a vertex, the corresponding time has to be integrated over, so we get a / dr, 
while for a closed loop we get an integral over the internal spatial momentum / (Pp/{2Tr)^. 
Considering the form of Gr we can already exclude the presence of loop with mixed lines, Hke in figure 2. Indeed, 




FIG. 2: The loop is formed by a retarded propagator Gr starting and ending at time rs. It is identically zero. 

such a loop would close a retarded propagator Gr on the same time t = ti = T2 but, due to the embedded 9{ti — T2), 
it vanishes. So the only possible loop that we can build with our set of Feynmann's rules is made of a full line. It is 
given by 

A(r) = I ^G^^\p,T,r), (22) 

where again by the superscript (0) we mean the free correlation functions. Since we will calculate correlation functions 
at higher orders, the superscript («) will help us to keep track of which order are we considering at each moment. As 
argued in [23|, primitively divergent graphs contain only vertices of the same type. If there were vertices of different 
type, then at least two internal lines would be retarded propagators, the corresponding momenta would be on shell, 
the corresponding loop integral would be finite and the graph would not have been primitively divergent. Now the 
graphs of the in-in effective action with all vertices of the same sign are just the graphs of the in-out theory plus their 
complex conjugates, so the primitive divergences must be the same. Once the primitive divergences are controlled, it 
is only a matter of combinatorics to show that the overlapping divergences disappear as well. 
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III. HIGHER ORDER CORRELATION FUNCTIONS FOR 



Since we are interested in the the IR modes, for which (— fcr <C 1), the free two-point functions Gc(k,Ti,T2) and 
Gnik, Ti,T2) can be expanded in powers of fcr (l^ : 

G(?)(fc,Ti,r2) = ^, (23) 



Gf{k,T,,T2) - e{T, - r2)^[k^{Tl - r|)]. (24) 

In the in-in formaHsm there are two vertices but we focus only on the one with three full lines and one dashed line. 
The reason is that the {4>c<t>c) = Gc has a momentum dependence fc~^ and thus is divergent in the infrared, while 
the {(j>c<l>A) = Gr does not. Moreover, we note that for a vertex with three dashed lines it is not possible to build 
loops since they vanish identically. 



A. First Order Diagrams 

The simplest correction to the free propagators Gq^ and G^^ comes from the graphs with a single tadpole. The only 
contribution comes from the graph with the full line loop thus figure [3] is the only first order correction to g[?^ : Due 




k~ ^ k- 



FIG. 3: One tadpole Gc and Gr propagators 

to the simmetry of Gc we must consider also the mirror diagram and using the expressions l(23l) and (|24l) , we obtain 
integrals of the form 

£[dT3a\T,)Gfik,n,T,) J ^G^")(p,T3,r3)G[?^(ft,T3,T2), (25) 

where p is the internal momentum of the tadpole and k the momentum flowing in the diagram. We set the inferior 
limit of integration to — instead of — oo, because we are interested in following perturbations from the moment of 
horizon exit up to some later time t. The horizon exit time is given by the condition — fcr/j ~ 1 and so r/j = — p The 
tadpole integral over d^p is divergent but can be in general regularized choosing appropriate infrared and ultraviolet 
cutoffs. Am and Ajjv, 



In our context the choice of the cut-offs is rather natural. The IR cut-off Aju is proportional to a^ff, where is the 
scale factor at the beginning of inflation and H is the Hubble rate, while Auv is equal to k, therefore the logarithm is 
proportional to the total number of e-folds from the beginning of inflation to the time the mode k exits the horizon. 
Before performing the calculation we must also consider the coefflcient in front of the graph coming from Wick's 
theorem. We have 0a(0c)"^ from the vertex and the external legs, (I)c{ti) and 4)c{t2)- So there are three possibilities 
for contracting 0(t2) with one of the (j)c of the vertex and one for contracting (j)c{Ti) with the vertex's (j)A, which 
sum up to 3 in front of the graph. Performing the calculation and considering also the mirror graph, at leading order 
we obtain 

Gg'(fc,Ti,r2) ^ ^^(ln(-fcri)+ln(-fcr2)). (27) 
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The retarded propagator at one loop has no mirror graph due to the oddness of the Gr under exchange of times. The 
Wick contraction factor is again 3. Then, at leading order for a one-loop Gr graph with ri > T2: 

-iG<^^\k,r,,r,) = 0ir,-r,)'^^in^ + r,^)\n^^^ . (28) 

These results coincide with those in Refs. and show the IR divergences due to ther cumulative effects. 



B. Higher Order Corrections 



To try to cure the divergences we need to proceed to higher orders. Already at the second order, three graph 
topologies can be identified, the tadpole chain (e.g. fig. S]), the tower graphs(e.g. fig. (H) and the sunrise (e.g. fig. [6]). 



1. Tadpole Cham Graphs 



The two-tadpole chain graphs can easily be calculated, basically adding a Gr and closing in a tadpole two of the 
straight lines of the second vertex (fig. The amplitudes for the two-tadpole chain propagators are 

G[?)(fc,n,r2) c ^(^"\ ^^{Xn\-kn)+\n\^kT2)), (29) 



The n-th order chain graphs can be calculated in the same way and we find 




FIG. 4: Two-tadpole Gc and Gr propagators. 



G^c?-"(fc,Ti,r2) 



2fc3 



(31) 
(32) 



- zG5j'^'"(fc, n, r2) = 0(ti - r2)^^(Ti3-V2^ - n^a^-^), (33) 

where e = ^ . The IR modes resummation of the chain diagrams generates a spectrum of perturbations which is no 
longer flat, but blue tilted. Next, we consider the resummation of the tower graphs. 



2. Tower Graphs 



To write the amplitude for the G^™'' at second order we must be more careful. Indeed, we have two loops 
which are chained one into the other; in the Gr tower diagram from two consecutive retarded propagators 
G^°Hfc,ri,T3)G(?)(fc , T3,T2) one obtains ti > T3 > T2 while for the Gc tower diagram one obtains only Ti > tj,. 
We note however that the times internal to the loops do not receive constraints from the 9 functions relative to 
Ti,T3,T2. The integral over the time dr4 must then extend from a loop characteristic time to the upper end which is 
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given by the 6 function embedded in Gii{p,T^,TA). The only time scale available is the one given by the momentum 
p, thus the integral over dT4 is evaluated between — i and T3. With these considerations the amplitude for the tower 
Gc diagram is given by 

{~tXa\Ts)) g[?^ {k, n,T2)H)Gf {p, T3, r4) 

{~tXa\n)) G^°\q,r4,n)G^°\p,T4,Ts), (34) 
and the amplitude for the retarded propagator is 

{-i\a\T:,)) Gf {k, T3, T2){-i)G^°\p, T3, r^) 

{-iXa\n)) G^^\q,n,n)G^°\p,n,Ts), (35) 
Performing the calculation and considering the mirror graph we obtain 

^towe.(2) ^ g (^)'i[ln^(-An)+W(-fcr,)], (36) 

while for the retarded propagator we obtain: 

Also in this case it is possible to calculate the rt-th order contribution: 

-zGr^("^(fc, n, r^) = ein ~ -2)^^'^" (^TH! (-1' + 1^'^" (|) ' 

and summing over all contribution with m between and 00 one obtains a contribution proportional to the one-loop 
tadpole diagram. For example, for the retarded propagator: 

-iG5?-'='-(fc,ri,r2) = (?(ri-T2)^(Ti3 + T23)4eln(^^^ . (38) 

Performing this last sums is equivalent to sum vertically over the whole class of tower graphs. Interestingly, the result 
is proportional to the first order graph, see eq. (28), that can be resummed as shown in the previous section. Therefore 
the resummation of the tower graphs can be accounted for just sending e into e' — be, which does not change the 
properties obtained from eq. l(3T|). 




FIG. 5: Second order tower diagrams for Gc{left) and GR{right). 
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FIG. 6: Second order sunrise diagrams for Gc{left) and GR{right). 
3. Sunrise Graphs 

Despite the graphical difference in respect to the tadpole chains, these diagrams translate exactly into tadpole graphs. 
For example the sunrise Gc can be written as 

{~iXa\T3)) G^^\p, T3, n){-i)Gf {k-p-q, T3, r4) 

{~iXa\T^)) G(?^(g,T3,r4)G[?^(fc,T4,r2), (39) 

and the same for the Gr. However the the combinatorial coefficient in front of the integrals is different than in the case 
of chain graphs, since there are two ways to contract 4>c(.ti) with a 4>/\ of the vertices, then three to contract one 4>c 
of the first vertex to the of the second and finally six to contract the remaining free 4>c in the two vertices. In total 
we have a 36/2!. The contributions coming from the tadpole and sunrise diagrams differ for a numerical constant. If 
we resum the whole class of sunrise diagrams the result is therefore proportional to the first order tadpole graph. One 
should not claim victory too soon though. Already at one loop, one should account for the vertex renormalization. 
Unfortunately, the IR resummation of the vertex renormalizing graphs proves to be a difficult task because of the 
presence of diagrams whose time fiow in the internal lines is not continous. This is not surprising at all, since it is 
very well known that only in certain class of self-interacting models, the full resummation is possible. We now turn 
to one of these examples, generalizing our results to iV ^ 1 fields respecting an 0{N) symmetry. 

IV. 0{N) SYMMETRY AND THE TRISPECTRUM 

One way to be able to discard all the graphs except towers and tadpoles is to assume that we have N :S> 1 fields 
with an 0{N) symmetry. Under this assumption, thanks to the normaHzation of the vertices (that gives a N^^ for 
each vertex) and the loop over the free indices in loops (an N for every free index) , tadpole and tower graphs are 
proportional to N°, while all the other graphs scale at least as N~^. Therefore, under the assumption that ^ 1 
and according to our diagrammatic results, we conclude that IR effects may be resummed. In Ref. [i^l the same 
model was analyzed using non perturbative stochastic techniques. There, from the Fokker-Planck equation, the field 
quadratic mean value, which is linked to the correlation functions by (0^) = G'+"'"(a;, a;),was obtained 

, , Tanh(^lna) , , 

--G++{x,x)^ -, (40) 

V A 




with A — 47r^A/3, and then inserted it into the gap equation in order to be able to solve for G++(a;,a:'), 

- + + ^ [cj^lix) + G++{x, x)] ) G++(x, x') = 



(41) 



The result for the two-point correlation functions was 



Gc{k;TuT2) « ^i~kT,)\-kT2)', (42) 



Gfl(fc;Ti,T2) « 0(ti - T2)^ [i-kTi)'{-kT2f-' ~ {- kT2)\- k^f ' ' )] . (43) 
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which have the same form as our eq. l(3T|) and l(33|) . except for the S = m^p/3i?^. The mass rrinp is what controls the 
IR divergences and appears because G+"'"(x, x) goes rapidly to a constant and thus plays the role of a mass in eq. (|40l) . 
In our diagrammatic approach, the same non perturbative mass appears if we think to eq. (22) as a gap equation for 
the IR cut-off, but using the resummed propagator G^^'", 



3AA : 



9XH^ 



that is 



rx 



2-K 



1 /A 



2 ' ' 2^ V 2 ' 



(44) 



(45) 



which has the same dependence on A as the non perturbative mass found in [2l||. The numerical coefficient do not 
coincide due to the different normaHzation of the potential. Indeed, redefining A A/3!, one obtains the potential 
V = A0^/4 and rup and e coincide with mnp and 5 of Ref. [2l|- Therefore, the diagrammatic approach reproduces the 
findings obtained using the stochastic approach. In particular, if we are interested in the evaluation of the trispectrum, 
we need to evaluate 



{5^c (ki , t)54>c (k2 , t)54>c (kg , T)5(pc (k4 , r)) = T(ki MM, ^i){2nf5 ^ (ki + ka + kg + k4) . (46) 
At tree level the 4-point function is built with free propagators as shown in figure ([7]) It is necessary to sum over the 




FIG. 7: Four point function calculated at tree level {left) and with the resummed propagators {right) 



permutations, since each of the four momenta can be fiowing through the retarded propagator. The amplitude at tree 
level is 



r--(ki,k„k3,k4) = Y.fdr'e-^^^t..^Or' i^~^K\r')'^ 



(47) 



{-^Gf {K r, r')) n G^°^ (k, , r, r')e-Et. ^ 



which, for /ct^ ^ 1 (i = 1, • • • , 4): 

T*-°(ki,k2,k3,k4,T) 



24ivnr=ifci 



(48) 



Equation (|48l) reproduces the result obtained by Bernardeau et al. in [25|. Under the hypothesis of 3> 1, the 
trispectrum is obtained the exactly resummed propagators, that is the 'double-line' propagators on the right of fig. [7l 
The amplitude is 



r=^-"(ki,k2,k3,k4) = E r dT'e-'^^t-^""^-' f-t^a\T') 

i—1 ^ 



(49) 



(-zGg'^'"(fc„ T, r')) n G^c "'"(fcj- , r, r')e-'Eti k,]r' 



X 1 



A^ 24 n"* k ^ ^ ' 

nA:/^^^ [-^4-4e(ifctT) + £;i_2e(*fctT)] 



(50) 
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where the function En{z) is defined as 

/•oo -zt 

En{z) = —dt. (51) 
Thus in the limit — fcr^ ^ 1 the amplitude becomes 

\ 77^4 y^4 ,3 

which again coincides in form with what found in Ref. j2l|. In particular, the resummed trispectrum shows an 
enhancement factor ~ l/VX compared to the tree level result. Higher loop corrections to the trispectrum are suppresed 
by the fact that the propagators are now IR regulated. 



V. CONCLUSIONS 



In this paper we have analyzed the IR corrections to the correlators for a self-interacting scalar field in a de Sitter 
background. We have used a full diagrammatic approach within the in-in formalism. In this sense, our results should 
be considered complementary to alternative approaches, e.g. the stochatic approach. It is reassuring that the same 
results are obtained once the IR resummation can be performed as in the large N 0{N) theory. 
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